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EXECUTIVE  SUMMARY 


The  quasi-static  antenna  design  algorithm  uses  multipole  basis  function  to  model  the  general  thick 
top  load.  A  sequence  of  solutions  converge  in  shape  and  Q.  The  absolute  minimum  Q-factor,  1.825, 
is  obtained  for  a  thick  disk  top  load  enclosed  by  a  sphere.  This  is  significantly  smaller  than  the  thin 
disk  top  load  Q-factor  2.349  and  previously  derived  thick-disk  Q-factor  2.078.  An  analytic  potential 
is  derived  for  each  multipole  basis  function.  The  capacitance  and  effective  height  is  calculated  from 
the  potentials  on  the  enclosing  sphere.  The  impedance  is  computed  with  Computer  Simulation 
Technology  (CST)  Microwave  Studio.  The  impedance  data  is  numerically  fit  to  a  dipole  eigenmode 
equivalent  circuit.  The  radiation  resistance  does  not  fit  the  expected  of  frequency  dependence 
(effective  height). 

The  error  is  an  of  term  that  is  explained  by  a  capacitor  approximated  for  the  octupole  eigenmode 
equivalent  circuit.  The  quasi-static  antenna  design  algorithm  predicts  the  DC  capacitance  and  the 
dipole  eigenmode  effective  height.  The  octupole  eigenmode  increases  the  radiation  resistance.  The 
Q-factor,  1.77,  is  lower  than  expected.  These  results  are  compared  to  the  spherical  cap  top  load.  The 
existence  of  an  of  limits  the  accuracy  of  the  theoretical  limits  in  the  Q  values  for  antennas. 


CONTENTS 

1.  INTRODUCTION . 1 

2.  CAPACITANCE  FOR  THE  OBLATE  SPHEROIDAL  ANTENNA . 5 

3.  GENERAL  CALCULATION  OF  RADIATION  RESISTANCE . 9 

4.  GENERAL  OBLATE  SPHEROIDAL  ANTENNA . 1 1 

5.  SECTION  CST  MODEL  OF  OBLATE  SPHERICAL  ANTENNA . 13 

6.  CST  MODEL  OF  A  SPHERICAL  TOP  LOAD . 19 

7.  CONCLUSION . 23 

8.  REFERENCES . 25 

Figures 

1.  First  four  top-load  radial  multipole  basis  function . 7 

2.  A  sequence  of  antenna  designs  with  one  to  five  load  multipoles . 1 1 

3.  The  antenna  shape  near  the  enclosing  surface . 12 

4.  Electric  field  calculated  as  a  function  of  distance  from  the  antenna  top . 12 

5.  The  equivalent  circuit  for  the  first  and  second  eigenmode  approximation . 13 

6.  CST  data  and  equivalent  circuit  model . 15 

7.  Error  in  eigenmode  model.  Typical  error  is  a  few  milliohms . 15 

8.  ft)2  and  co4  contribution  to  resistance . 16 

9.  Q  for  the  minimum  Q  OSA  design . 16 

1 0.  Q-factor  for  the  minimum  Q  OSA  design . 17 

1 1 .  CST  impedance  and  two  eigenmode  model . 20 

1 2.  Error  in  the  two-eigenmode  model . 21 

13.  ft)2  and  co4  components  of  resistance . 21 

1 4.  Q  for  spherical  top  load . 22 

15.  Q-factor  for  spherical  top  load . 22 

Tables 

1 .  Least  square  best  fit  to  Stuart's  eigenmode  impedance  circuit  model . 14 

2.  Spherical  top-load  data  fit  to  a  two-eigenmode  model . 20 


v 


1.  INTRODUCTION 


Electrically  small  antennas  are  common  to  portable  electronics  devices,  personal  digital  assistants 
(PDAs),  cell  phones,  etc.  These  antennas  are  small  compared  to  wavelength  or  electrical  small 
antennas.  Designing  electrically  small  antennas  is  an  art;  they  normally  have  a  low  radiation 
resistance,  a  large  reactance,  and  a  small  bandwidth.  Additional  matching  components  are  required  to 
eliminate  the  reactance  and  increase  the  resistance.  The  objective  of  this  work  is  to  increase  the 
radiation  resistance,  decrease  reactance,  and  maximize  the  bandwidth.  The  quasi-static  antenna 
design  algorithm  is  more  science  than  art. 

L.  J.  Chu  [1]  derived  a  lower  limit  for  the  Q,  for  electrically  small  antennas.  Chu's  Q  calculation 
is  based  on  the  radiated  energy  and  the  stored  energy  outside  a  sphere  enclosing  the  antenna;  the 
energy  inside  the  sphere  is  assumed  to  be  zero.  H.  L.  Thai  [2]  refined  this  limit  by  assuming  the 
antenna  current  is  limited  to  the  surface  of  the  enclosing  surface.  The  folded  spherical  helix  design  by 
S.  R.  Best  [3]  meets  Thai's  limit  for  TM10.  M.  Gustafsson,  C.  Sohl,  and  G.  Kristensson  [4]  derived  a  Q 
limit  based  on  the  optical  theorem.  A.  D.  Yaghjian  and  H.  Stuart  [5]  derived  a  more  restrictive  limit  on  Q . 

The  energy  inside  the  sphere  limits  the  antenna  performance.  The  thin  disk-loaded  monopole  is  a 
leaky  capacitor  with  a  large  amount  of  stored  energy  inside  the  sphere.  The  electric  field  under  the 
disk  is  larger  than  the  field  above  the  disk.  The  quasi-static  antenna  design  algorithm  [6,  7] 
analytically  computed  a  thick  disk-loaded  dipole.  The  thick  disk-loaded  dipole  reduces  stored  energy 
and  Q  by  filling  some  of  the  spherical  volume  with  conductor.  The  electric  field  under  the  top  load  is 
still  larger  than  the  field  above  the  top  load.  This  paper  shows  how  to  reduce  the  electric  field  under 
the  top  load  and  increases  the  radiation  resistance  of  the  antenna  [8],  The  antenna  shape  shifts  to 
eliminate  the  gap  (and  the  energy)  between  the  top  of  the  antenna  and  the  enclosing  sphere.  The 
stored  energy  inside  the  enclosing  sphere  is  limited  to  the  region  between  the  antenna  and  the  ground. 

Electrically  small  antennas  have  electric  fields  much  larger  than  the  magnetic  fields  below  the 
antenna  resonance.  The  Asymptotic  Conical  Dipole  (ACD)  [9,  10,  11]  was  the  first  antenna  design 
with  quasi-static  methods.  In  electrostatics,  a  perfect  conductor  is  the  same  as  an  equipotential 
surface.  A  line  of  constant  charge  on  the  z-axis,  with  an  image,  will  generate  the  ACD  antenna 
design.  Each  ACD  antenna  has  a  different  height.  The  quasi-static  antenna  design  algorithm  [7]  fixes 
the  antenna  height  to  a  constant  a  and  the  length  of  the  line  charge  m  is  varied;  the  antenna  fits 
within  an  enclosing  sphere  with  a  radius  a.  The  parameter  k  is  a  dimensionless.  The  Q  is 
calculated  from 


Q  = 


l 

ojCRKad 


(1) 


where  RrcuI  is  radiation  resistance,  C  is  capacitance,  and  co  is  angular  frequency.  The  quantities 
RRad  and  C  are  functions  of  k.  The  radiation  resistance  is  calculated  from  the  effective  height.  The 

capacitance  is  calculated  from  the  charge  on  the  antenna  arm  and  the  maximum  potential  on  the 
spherical  enclosing  surface  with  a  radius  a. 
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The  above  equation  for  Q  is  valid  below  resonance  and  it  gives  only  the  first  term  in  Chu's 
equation: 

Q>Jkaf  +  Ja,  (2) 

where  k  =  2n  t  X  and  a  is  the  radius  of  the  enclosing  sphere.  The  ACD  design  is  extended  by  adding 
a  disk-shaped  charge  distribution  as  a  load  on  the  line  charge  [7],  A  disk  in  free  space  is  used  as  the 
charge  distribution.  The  disk  charge  distribution  is  symmetric.  The  electric  fields  between  the  disks 
and  image  will  be  larger  than  the  electric  fields  above  the  disk.  There  is  no  requirement  for  the  disk 
charge  distribution  to  be  symmetric.  Adding  a  dipole  moment  moves  the  charge  from  the  bottom  of 
the  disk  to  the  top  of  the  disk.  This  reduces  the  electric  field  between  the  disks.  A  series  of  multipole 
charge  distributions  can  be  added  to  the  disk  to  model  the  general  charge  distribution  on  the  disk. 

Sections  2  and  3  are  not  needed  to  understand  the  results  given  in  Sections  4  and  5.  Section  2 
shows  how  the  potential  and  top  load  capacitance  is  computed  from  electrostatic  solutions  in  oblate 
spheroidal  coordinates.  Each  solution  represents  a  unique  multipole  moment  with  an  unique 
potential.  Only  the  rotationally  symmetric  top-load  multipole  modes  will  be  used  in  this  model.  The 

n ,h  multipole  falls  off  as  in  the  far  field.  The  potential  for  the  dipole  multipole  term  and  high 

odd  moments  add  in  the  far  field.  The  monopole  and  higher  even  moments  cancel  in  the  far  fields. 

All  of  the  multipole  make  a  unique  and  diminishing  contribution  to  the  near  fields. 

The  perfect  conductor  boundary  condition,  =  0,  requires  the  charge  distribution  to  be  enclosed 

by  the  antenna  surface.  Only  a  subset  of  the  charge  distributions  satisfies  this  boundary  condition.  The 
multipole  moments  have  negative  potentials,  which  can  cause  the  equipotential  surface  to  terminate 
on  the  disk  or  feed  wire.  This  requires  an  additional  step  in  the  solution  process;  the  equipotential 
surface  is  sampled  to  verify  that  the  charge  is  enclosed  by  the  equipotential  surface.  The  final 
solution  must  be  verified  with  a  detailed  calculation  of  the  antenna  shape.  In  Section  3,  the  effective 
height  is  calculated.  The  effective  height  calculation  with  dipole  moment  top  load  is  non-trivial.  The 
effective  height  cannot  be  calculated  with  the  conventional  formula  [8],  The  effective  height  is 
calculated  indirectly  from  the  potential  on  the  enclosing  sphere. 

In  Section  4,  a  sequence  of  multipole  basis  functions  are  used  to  design  minimum  Q  antennas.  The 
antenna  designs  appear  to  converge  in  both  shape  and  Q.  The  final  antenna  design  fills  the  top  of  the 
sphere;  the  area  under  the  antenna  is  the  only  region  with  electric  fields  that  contributes  to  stored 
energy  inside  the  enclosing  sphere.  The  electric  field  on  the  surface  of  the  antenna  is  plotted.  The 
electric  field  under  the  antenna  is  reduced  by  adding  multipole  moments.  The  final  antenna  also  has 
an  almost  horizontal  lower  surface. 

In  Section  5,  Computer  Simulation  Technology  (CST)  Microwave  Studio  is  used  to  calculate  the 
impedance  and  Q  for  a  1-m  high  antenna  with  a  1-cm-diameter  hollow  feed  line.  The  impedance  is 
calculated  with  a  sequence  of  energy-based  adaptive  iterations.  L.  J.  Chu  [1];  H.  D.  Foltz,  J.  S. 
McLean,  and  L.  Bodner  [12];  and  H.  R.  Stuart  [13]  used  a  series  L ,  C  with  a  resistor,  R0 ,  parallel 

to  the  inductor  to  model  both  the  reactance  and  co2  dependence  radiation  resistance.  The  model 
represents  the  dipole  eigenmode.  The  least  squares  fit  to  the  CST  impedance  is  a  good  fit  to  the 
reactance.  The  radiation  resistance  has  an  obvious  error.  The  radiation  resistance  has  an  co1  and  co4 
frequency-dependence  terms;  the  coefficients  are  calculated  with  a  least  square  fit.  The  co2 
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contribution  gives  the  same  effective  height  as  the  quasi-static  antenna  design  algorithm.  The  co4 
frequency-dependence  term  is  the  source  of  the  error  in  the  circuit  model. 

The  next  eigenmode,  the  octupole,  is  a  circuit  connected  in  parallel  with  the  first  eigenmode  model 
[13].  The  octupole  (and  higher)  equivalent  circuits  reduces  to  a  capacitor  at  low  frequencies.  This 
capacitance  combined  with  dipole  eigenmode  circuit  gives  an  anti-resonance.  The  dipole  and  octupole 
eigenmode  currents  interfere  with  each  other  to  create  a  high  impedance  at  anti-resonance.  The 
anti-resonance  does  not  have  a  deep  physical  meaning.  The  combined  dipole  eigenmode  model  and 
higher  eigenmode  capacitance  significantly  improve  the  numerical  fit  to  the  CST  impedance  data. 

The  capacitance  value  was  calculated  from  the  anti-resonance  frequency;  the  other  circuit  elements 

L,  C,  and  -^o  ,  were  unchanged.  The  improved  fit  to  the  radiation  resistance  data  is  attributed  to  a 
co 4  frequency-dependant  term  introduced  by  the  parallel  capacitance.  The  simple  equivalent  circuit 
for  the  octupole  does  include  any  radiation  resistance.  The  higher  radiation  resistance  reduces  the  Q  - 
factor  below  the  value  predicted  by  the  quasi-static  antenna  design  algorithm.  The  extra  capacitance 
introduced  by  the  feed  line  is  estimated.  The  quasi-static  antenna  design  algorithm  gives  both  the  DC 
capacitance  and  the  co1  dipole  radiation  resistance. 

In  Section  6,  CST  is  also  used  to  calculate  the  impedance  of  the  spherical  top-loaded  monopole. 

The  same  feed  wire  was  used  for  the  1-m  spherical  top-loaded  monopole.  The  data  is  fit  to  the  same 
dipole  eigenmode  with  a  capacitor  approximation  of  higher  eigenmodes.  The  capacitance  of  the  top¬ 
load  is  larger  than  the  above  design.  The  effective  height  is  smaller.  The  Q-factor  is  lower  than  the 
result  obtained  by  A.  R.  Lopez  [14].  Electrostatic  data  is  not  available  to  compare  to  the  circuit  results. 

Section  7  is  the  conclusion.  The  quasi-static  antenna  design  algorithm  gives  a  solution  very  close 
to  the  spherical  capped  monopole  (the  top-loaded  design  with  the  lowest  Q).  The  algorithm  does  not 
model  all  possible  antenna  shapes.  The  algorithm  does  not  allow  the  bottom  surface  of  the  antenna  to 
be  above  the  source  charge  on  the  disk.  A  thick  spherical  cap  top-load  antenna  could  have  a  lower  Q. 
The  quasi-static  antenna  design  algorithm  produces  a  very  good  antenna  design  with  modest  effort. 
Other  designs  can  be  computed  with  a  significantly  higher  radiation  resistance  and  a  slightly  larger 
Q.  The  algorithm  can  be  adapted  to  other  enclosing  surfaces. 
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2.  CAPACITANCE  FOR  THE  OBLATE  SPHEROIDAL  ANTENNA 


The  detailed  discussion  of  the  algorithm  is  given  in  Reference  [7] .  This  discussion  is  the  more 
general  case.  The  capacitance  of  the  general  thick  top  load  is  calculated  from  the  potential  on 
enclosing  surface  ( q  =  CV ) .  The  potential  is  calculated  from  an  analytic  solution  to  the  electrostatic 
equations.  The  top  of  the  hollow  feed  line  has  a  height  m  and  {\-r)q  charge;  the  potential  [15,  16] 
is 


tikh/MMl 

4  nejea  +  £,-)/ 


(3) 


where  Sm  =  is  for  the  monopole,  dj  =  is  for  the  image  monopole,  Rt  =  -J( z-m )2  +  p1  is 
the  distance  from  p ,  z  to  the  top  of  the  monopole,  r  =  ^ z 2  +  p 2  is  the  distance  from  p,z  to  the 
feed  point,  and  r  =  J(z+mf  +p2  is  the  distance  from  p,  z  to  the  bottom  of  the  image  monopole.  The 
parameter  r  is  the  fraction  of  the  total  charge  q  on  the  antenna  top  load.  The  superscript  indicates 
an  ACD  stem.  The  parameter  k  is  dimensionless  and  0  <  k  <  1 . 

The  disk  also  has  the  height  m  ,  a  disk  radius  of  aa,  and  a  net  charge  zq  where  the  parameter  a 

is  dimensionless  and  has  the  range  of  values  0  <  a  <  ^1-k1  .  G.  Arfken  [17,  pp.  599-601  and  596] 
gives  the  general  solution  to  potential  on  the  disk  V2d>(»,v,0)  =  0.  The  potential  on  the  disk  is  a 
linear  combination  of  multipole  moments  Kn(u)P"'( cosv)e""^  where  v  and  u  are  defined  by  the 
oblate  spheroidal  coordinate  system: 


z  =  aa  sinh(w)  cos(v), 

(4) 

x  =  oca  cosh(w)  sin(v)  cos(^), 

(5) 

y  =  aa  cosh(w)sin(v)sin(^), 

(6) 

where  the  aa  is  the  disk  radius..  The  antenna  problem  will  be  simplified  by  assuming  rotational 
symmetry  where  p  =  aacosh(w)cos(v)  : 

oo  m=n  i 

®(».V,#)=  ZZ v)e-R  (7) 

n= 0  m= 0  ^ 
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The  rotational  symmetry  in  (j)  requires  bmn  =  0  for  m  >  0  and  the  associated  Legendre 
Polynomials  reduce  to  Legendre  Polynomials,  Tf  (cosv)  =  Pn (cos  v) .  The  1  la"  is  added  for 
numerical  convenience.  The  Kn{u )  functions  are  as  follows: 


if0(sinhit)  =  arccot(sinhtt), 

(8) 

ifj(sinhw)  =  l-sinh(i/)arccot(sinh(w)), 

(9) 

K2  (sinh  u )  =  [(3  sinh 1  (u).  + 1)  arccot(sinh(«))  -  3  sinh(w)  J/  2, 

(10) 

and 

^5  3  3  5  63  j  5 

^(sinht/)  = - sinh3(it)  +  — sinh(it)  arccot(sinh(w))  +  —  — sinh2(w)  +  — 

y2  2  J  3y2  2j 

1 
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(ID 

The  general  case  is 

(12) 

where 

QM)-(  l)Pn('Z)*K  COt(^)  ^  L  }Pn-Al%)  - 

77  3(77  -  1 ) 

(13) 

The  function  Kn  at  large  distance  is  calculated  by  expanding  a  cot(snh(u))  =  a  tan(l  /  sinh(w))  in  a 


power  series  of  - - .  After  simplifying,  the  leading  term  is 


sinh(w) 


A'o  fts 


1 


sinh(u) 


for  u»  2. 


(14) 


and 


K, 


2  1 
3  sinh2(w) 


for  u  »  2, 


(15) 


2  1 
15  sinh3(tt) 


for  u  »  2, 


(16) 


2  1 
35  sinh4(tt) 


for  it  »  2 


(17) 
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since  aasinh(w)  =  z ,  the  above  equations  reduce  to  — ,  i(^)2,  and  ^■(J?L)4-  This  is  the 

expected  multipole  expansion  of  the  field  on  the  z-axis.  Figure  1  plots  Kn  as  a  function  of  sinh(w); 
this  shows  the  dependence.  The  factor  —  is  included  in  Equation  7  to  eliminate  the  a" 

dependence  in  the  far  field.  The  value  of  b00  =  jj—  is  calculated  from  the  solution  of  the  disk  in  free 
space  V(u )  =  4Jaa  arccot(sinh  u) ,  where  rq  is  the  net  charge  on  the  disk.  The  only  restrictions 
placed  on  the  values  of  bnm  is  the  charge  distribution  must  be  enclosed  within  an  equipotential 
surface. 


Radial  Dependence  of  the  Disk  Basis  Functions 


sinh(u) 

Figure  1.  First  four  top-load  radial  multipole  basis  function. 

The  combined  potential  is 

0(z,  p)  =  cp  "CD  (p,  z)  +  O  TopDlsk  (ut ,  v( )  -  O  Rotmk  (uh ,  v6 ),  (18) 

where  the  variables  ut ,  v  t ,  and  u  t ,  v  h  are  defined  from 


z  —  kci  =  oca  sinh(w; )  cos(v, ), 

(19) 

p  =  «acosh(w,)sin(v,), 

(20) 

z  +  Ka  =  aa  sinh(wA )  cos(vA ), 

(21) 

p  =  «acosh(wi>,)sin(vz,). 

(22) 
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The  capacitance  is  calculated  from  the  maximum  value  of  ® {z,p)  on  the  sphere  z  =  a  cos  6  and 
p  =  a  sin  0 : 

=  max(®(a  cos  0 ,  n  sin  0  j  j  where  0  ^  0  ^  7t  /  2  (23) 


and 


C  = 


® 


(24) 


Max 


In  the  above  coordinate  systems,  vt  =  0  and  vb  =  0  are  defined  in  the  +  z  direction.  The  charge 
distribution  that  generates  the  term  0Bo,Dlsk  (iib,vh)  is  the  mirror  image  of  the  charge  distribution  for 

®TopDisk{u„vt): 

00  h 

<bT“’Di“(u„v,)  =  '£^K„(u,K(c  osv,)  (25) 

n= 0  a 


and 


® ' ““  K  V, J  =  Z  (- 1 )"  h K„  (“» « (cos  V,  )  . 


(26) 


The  sign  on  prKn (ub)Pn (cos  )  follows  from  Pn  (cos  vb ) ,  being  even  or  odd  relative  to  v  =  n  /  2 . 

The  dipole  term  is  the  simplest  example  to  understand:  Kt (ut)Pt (cos  vt)  has  positive  charge  on  the 

top  of  the  disk  and  Kt(uh)Pl(cosvh)  has  positive  charge  on  the  top  of  the  image  disk  (top  is  +z 

direction  for  both  disks).  Both  terms  have  a  negative  charge  on  the  bottom  of  the  disk  and  image 
disk.  The  combined  terms 

Kx  (ut  )PX  (cos  v,)+Kx  (ub  )PX  (cos  vb )  (27) 


are  odd  about  z  =  0.  For  the  case  bx  >  0 ,  the  dipole  moment  of  the  antenna  is  increased  and  the 
electric  field  under  the  antenna  is  decreased  (with  the  other  parameters  unchanged).  In  the  general 
case,  when  n  is  odd,  the  terms 

Kn(ut)Pn(cosvt)  +  Kn(ub)Pn(^svb)  for  odd  n  (28) 

increase  the  n,h  moment.  When  n  is  even,  the  terms 

K„  (u,  )P„  (cos  v, )  -  Kn  (uh )Pn  (cos  v, )  for  even  n  (29) 

subtract  in  the  far  field  the  first  contributing  far-field  moment  is  the  (n  + 1  )'h  moment.  Each  term  has 
a  unique  contribution  to  the  electric  field  under  the  antenna. 

The  next  step  is  to  compute  the  radiation  resistance. 


3.  GENERAL  CALCULATION  OF  RADIATION  RESISTANCE 


The  effective  height  for  a  rotationally  symmetric  charge  distribution  is 

1  pz=a  rp=^J a 2  — z 


h  = 

nEff 


i  rz=“  r/ 

qNe,  •’-"=0  ^=0 


zq{p,  z)2npdpdz. 


(30) 


From  a  previous  paper  [6,  7],  this  simplifies  to 

m 


hEfr=—(\-T)  + 


/car. 


(31) 


How  to  apply  the  above  effective  height  calculation  to  the  dipole  term  K]  (a)P]  (cos  v)  is  not 
obvious.  The  ^j(m)TJ(cosv)  term  is  discontinuous  at  u  =  0.  The  function  is  odd  about  v  =  k  /  2  or 
Kx  (O)T;  (cos  v)  =  -Kx  (0)7)  (cos(;r  -  v)) .  The  dipole  term  would  give  zero  for  a  thin  disk  with  no  net 
charge.  The  above  equation  is  only  valid  for  even  potentials  on  the  disk  (monopoles,  quadrupoles, 
etc.). 

To  derive  the  general  case,  the  dipole  moment  used  as  the  starting  point  is 


/  2  2 
rz=a  rp=-y  a  -z 

=  zq(p,z)2npdpdz. 

Jz=-a  Jp=0 


(32) 


The  potential  outside  (and  on)  the  sphere  can  be  represented  as  a  sum  of  spherical  harmonics: 


a ^ 

J-) 


n+ 1 


7)(cos£)  for  r>< 


(33) 


where  pz  =  fiwa  and  r  =  ^ z 2  +  p2  with  all  of  the  charge  enclosed  in  a  sphere  of  radius  a  .  The 
dipole  moment  is  computed  from  the  potentials, 

3 


p.  =  —a 
2 


£  q  ®(a,6,)T|(cos6,)sin$7#. 


(34) 


This  term  is  independent  of  the  charge  distribution.  The  effective  height  is 


hEff  ~ 


1 


2  q 


P_ 


(35) 


Net 


9 


4.  GENERAL  OBLATE  SPHEROIDAL  ANTENNA 

In  the  calculation  for  the  thin  disk-loaded  monopole,  the  feed  line  is  neglected.  The  capacitance 
calculation  assumes  all  of  the  charge  is  on  the  top  load  disk.  This  model  includes  a  hollow  feed  line. 
The  Matlab B  calculation  of  the  Oblate  Spheroidal  Antenna  (OSA)  assumes  that  99%  of  the  charge  is 
on  the  top  load  and  a  token  1%  of  the  charge  is  on  the  feed  line:  z  =  0.99  and  (l  -z)  =  0.01 .  As 
shown  in  Figure  2,  the  antenna  designs  appear  to  converge  in  shape  and  Q . 

Antenna  Cross  Section  for  Different  Moments  with  x  =  0.99 


Figure  2.  A  sequence  of  antenna  designs  with  one  to  five  load  multipoles. 

Figure  3  shows  a  detailed  plot  of  the  antenna  surface  near  the  sphere.  Spherical  coordinates 
r  =  Jz2  +  p1  are  used  for  the  vertical  axis  and  0  =  arcsin {pi r)  is  used  for  the  horizontal  axis.  As 

basis  functions  are  added  the  antenna  design,  the  shape  approaches  the  surface  of  the  enclosing 
sphere.  Figure  3  shows  that  adding  extra  multipole  moments  has  a  diminishing  impact  on  the 
Q-factor.  The  multipole  moments  fall  off  as  l/rn+l ;  they  have  diminishing  impact  on  the  electric  field 
outside  the  sphere  and  under  the  antenna.  The  electric  field  is  calculated  from  the  potential  and 
plotted  in  Figure  4. 

The  vertical  scale  is  electric  field  magnitude  plotted  as  a  function  of  distance  along  the  surface 
measured  from  the  top  of  the  antenna,  z  =  1  and  p  =  0  .  The  spike  in  the  electric  field  is  at  the 
transition  from  almost  spherical  surface  to  the  underside  of  the  antenna.  For  the  solution  with  five 
multipoles,  this  transition  is  very  sharp.  The  electric  field  under  the  antenna  is  reduced  by  the  dipole 
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moment  term.  Figure  4  shows  that  the  position  of  the  spike  moves  to  the  right  as  multipoles  are 
added.  This  increases  the  surface  area  of  the  top  surface  of  the  antenna  and  the  total  charge  on  the  top 
of  the  antenna. 


Antenna  Cross  Section  Near  Inclosing  Sphere  with  z  =  0.99 


Figure  3.  The  antenna  shape  near  the  enclosing  surface. 


Electric  Field  as  a  Function  of  Distance  z  =  0.99 


Figure  4.  Electric  field  calculated  as  a  function  of  distance  from  the  antenna  top. 
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5.  SECTION  CST  MODEL  OF  OBLATE  SPHERICAL  ANTENNA 

Computer  Simulation  Technology  (CST)  Microwave  Studio  was  used  to  model  the  antenna.  An 
accurate  model  for  the  energy  is  critical  for  Q  calculations.  CST  includes  an  energy-based  adaptive 
meshing  that  increases  the  mesh  resolution  in  areas  of  high-energy  density.  A  sequence  of  adaptive 
meshing  models  were  run  to  refine  the  starting  mesh  for  the  next  model1.  A  4.3-ohm  cylindrically 
edge  port  was  used  at  the  hollow  feed  point  and  the  excitation  frequency  range  was  from  15  to  25  MHz; 
this  reduces  the  energy  reflected  at  the  input  to  the  antenna2 *.  The  hollow  feed  line  was  modeled  independ¬ 
ently  to  determine  the  required  meshing.  The  convergence  in  the  solutions  is  evaluated  by  looking 
at  the  largest  difference  between  previous  and  current  S  n  .  In  the  final  run,  this  difference  is  3.4e-4. 

The  process  was  time  consuming,  but  it  eliminated  the  risk  of  using  too  fine  a  mesh  and  introducing 
numerical  errors  in  the  solution. 

The  antenna  capacitance  and  effective  height  can  be  extracted  from  the  CST  impedance  data 
( p  <  0.5  )  with  a  MatLab  program.  The  MatLab  program  calculated  the  least  square  best  fit  to 
Stuart's  eigenmode  impedance  circuit  model,  Figure  5  and  Table  1.  The  fit  to  the  reactance  was  very 
good  but  the  radiation  resistance  has  a  large  error.  An  additional  octupole  eigenmode  can  be  modeled 
at  low  frequencies  as  a  parallel  capacitor  COSAOP+  where  the  subscript  indicates  octupole  plus  higher 
order  eigenmodes  (the  inductors  short  the  other  circuit  elements  at  low  frequencies).  This  extra 
capacitance  adds  an  antiresonance  to  the  single  eigenmode  impedance.  Another  CST  model  was  run 
to  calculate  the  antiresonance  frequency,  90.916  MHz.  The  value  of  COSAOP+  was  calculated  based  on 

this  frequency  and  the  previously  calculated  values  of  R0 ,  C ,  and  L  . 


Figure  5.  The  equivalent  circuit  for  the  first  and  second  eigenmode  approximation. 


Figure  6  shows  the  CST  impedance  and  the  two  eigenmodes  circuit  approximation.  The  difference 
is  so  small  it  is  plotted  in  Figure  7.  The  circuit  parameters  are  in  Table  1.  The  C^SJ  capacitance  is 


^ CST  _  ^  CST  (-i 

DC  ~  ^OSA  +  L'OSAOP+ 


95.032. 


(36) 


'Shrikrishna  Hegde  at  Sonnent  Software  provided  insight  into  energy  adaptation  and  meshing. 

2James  Whillhite  at  Sonnent  Software  suggested  matching  the  a  source  resistance  to  the  antenna 

and  using  a  narrow  frequency  range  to  maximize  the  energy  delivered  to  the  antenna. 
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Table  1.  Least  square  best  fit  to  Stuart's  eigenmode  impedance  circuit  model. 


Parameters 

Eiqenmode  1 

Eiqenmode  1  &  2 

QSADA 

*3  OSA 

90.2683  pF 

90.2683  pF 

90.364  pF 

OSA 

0.6758598  pH 

0.6758598  pH 

N/A 

^■OOSA 

1.96632  KQ 

1.96632  KQ 

N/A 

' RES 

20.39706  MHz 

20.39705  MHz 

N/A 

^Eff 

N/A 

0.7063 

0.7108 

*3  OSA  OP+ 

N/A 

4.7638  pF 

N/A 

Antiresonance 

N/A 

90.916  MHz 

N/A 

Q-factor 

N/A 

1.77 

1.82 

The  impact  of  the  finite  diameter  feed  line  can  be  estimated.  The  current  on  the  feed  line  is  almost 
constant;  almost  all  of  the  charge  flows  onto  the  top  load.  A  reasonable  approximation  for  the  charge 
distribution  on  the  feed  line  is  a  triangular  distribution.  The  net  charge  on  the  feed  line  is  1/2  the 
monopole  charge  (constant  charge  distribution  for  electrical  small  monopole).  The  capacitance  of  a 
monopole  0.587-m  high  is  6.2  pF.  The  feed-line  capacitance  is  about  3.1  pF.  The  quasi-static  antenna 
design  algorithm  (QS  ADA)  assumes  1  %  of  the  total  charge  is  on  the  feed  line  or  the  feed  line 
capacitance  is  0.9  pF  (=  0.01  *  90.364  pF) .  This  difference  is  CStem  =2.2  pF.  The  difference 
between  the  two  calculations  is  about  2.4  pF  or  about  2.6%: 

cg  =  cg,  +CS""=92.6.  (37) 

The  resistance  has  two  components,  of  and  co4  terms 

Rosa  =  bi (®  / ar  f  +  b2(a>  l  co,. )4 ,  (38) 


where  cor  is  the  angular  resonant  frequency,  bx  =  3.6465  and  b2  =  0.4575.  Figure  8  shows  the 

curve  fit.  The  effective  height  in  Table  1  is  calculated  from  the  of  term  in  the  radiation  resistance;  it 
agrees  with  the  quasi-static  antenna  design  algorithm.  The  radiation  resistance  circuit  element  is  not 
included  in  the  octupole  eigenmode  circuit  model.  The  co4  term  is  caused  by  a  circulating  current  in 
the  parallel  CL  circuit.  This  current  is  higher  than  the  current  in  the  dipole  eigenmode  equivalent 
circuit.  The  antenna  Q  will  be  lower  than  predicted  by  the  quasi-static  antenna  design  algorithm. 


The  Q,  Figure  9,  was  calculated  using  the  Yaghijian  and  Best  [18]  result: 


g  =  — 1 
2 


\{  dR  ^ 

2 

( 

x\ 

dX\ 

+ 

-  H - 

\doo  ) 

V 

(O 

dco  ) 

(39) 
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Figure  10  shows  the  Q(kaf  “Q-factor”  for  the  antenna  that  was  calculated  from  a  two-eigenmode 
circuit  approximation.  The  Q  curves  are  calculated  numerically  from  the  fit  data.  The  minimum 
Q-factor  is  1.77. 
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Figure  6.  CST  data  and  equivalent  circuit  model. 


Figure  7.  Error  in  eigenmode  model.  Typical  error  is  a  few  milliohms. 
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The  Dipole  and  Octipole  Contribution  to  Resistance 


Figure  8.  co1  and  co4  contribution  to  resistance. 


The  Q  for  the  Minimum  Q  Quasistatic  Design 


Figure  9.  Q  for  the  minimum  Q  OSA  design. 
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The  O  factor  for  the  Minimum  0  Ouasistatic  Design 


Figure  1 0.  Q-factor  for  the  minimum  Q  OSA  design. 
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6.  CST  MODEL  OF  A  SPHERICAL  TOP  LOAD 

The  same  method  was  used  to  model  a  spherical  top  load.  The  spherical  top  load  is  a  much  harder 
problem.  The  edge  of  the  spherical  top  load  will  have  a  much  higher  field  than  the  above  solution; 
this  requires  a  higher  resolution  near  the  edge  of  the  spherical  top  load.  If  the  edge  thickness  was  set 
too  small,  CST  will  not  accurately  model  the  field.  The  shell  thickness  is  0.83-cm;  the  horizontal 
edge  is  1  cm  and  the  height  of  the  edge  is  0.560  m.  (The  exact  minimum  is  at  0.556.)  The  feed  wire  is 
the  same  as  in  the  above  CST  model.  The  same  eigenmode  circuit  model  was  used  to  fit  the 
impedance.  The  circuit  parameters  are  listed  in  Table  2.  In  this  case,  the  capacitance  and  effective 
height  are  not  known.  Alfred  R.  Lopez  [14]  published  a  Q-factor  of  Q(kaf  -  1.75. 

The  same  Matlab  program  was  used  to  numerically  fit  the  data  to  a  two-eigenmode  model,  Table  2 
and  Figure  11.  Figure  12  shows  the  magnitude  of  the  error  is  <  0.0013  Q  .  If  the  least  squares  fit  was 
expanded  to  include  the  variable  Cop+ ,  the  error  would  be  smaller.  The  spherical  top  load  has  a 

larger  surface  area  and  a  smaller  charge  density  larger  than  the  above  antenna.  The  potential  on  the 
spherical  top  load  will  be  smaller.  The  DC  capacitance  is  8%  larger  than  the  above  antenna.  The 
effective  height  is  5%  smaller  than  the  above  antenna.  A  thinner  spherical  shell  will  have  less  charge 
on  the  edge  and  a  higher  charge  density  on  the  sphere.  The  higher  charge  density  increases  the 
potential  on  the  sphere;  this  decreases  the  capacitance.  Moving  the  edge  charge  to  a  higher  position 
increases  the  effective  height.  Stuart  [19]  pointed  out  that  a  thicker  shell  has  a  higher  Q-factor.  The 
feed-line  analysis  is  the  same  as  above.  A  1-m  monopole  has  9.54-pF  capacitance.  The  feed  line 
would  have  a  capacitance  around  5  pF.  The  CST  model  is  reasonable  accurate. 

The  radiation  resistance  has  two  components,  co1  and  a>A : 

RSP  =sx(gjI  (or  )2  +  s2  (co  /  cor  )4 ,  (40) 

where  cor  is  the  angular  resonant  frequency,  Sj  =  1.960684  and  s2  =  0.28059356.  Figure  13  plots 

the  total  radiation  resistance,  the  co2  and  co4  terms.  The  Q  is  plotted  in  Figure  14  and  the  Q-factor  is 
plotted  in  Figure  15. 
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Table  2.  Spherical  top-load  data  fit  to  a  two-eigenmode  model. 


Parameters 

Eigenmode  1 

Eigenmode  1  &  2 

Csp 

97.20870562  pF 

97.20870562  pF 

Lsp 

1.0303518  pH 

1.0303518  pH 

^0  SP 

4.7281347  K D. 

4.7281347  KD 

f  SP  RES 

15.90688  MHz 

15.90654  MHz 

hEff 

0.6649 

0.6649 

Cop 

N/A 

6.4012  pF 

f An  tires 

N/A 

63.964  MHz 

Q  -  Factor 

N/A 

1.70 

CO 

E 

-C 

o 


The  Spherical  Cap:  CST  Impedance  and  Two  Eigenmode 
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Figure  1 1 .  CST  impedance  and  two  eigenmode  model. 
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Ohms 


Figure  12.  Error  in  the  two-eigenmode  model. 


The  Dipole  and  Octipole  Contribution  to  Resistance 
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The  O  for  Spherical  Top  Load 


Figure  14.  Q  for  spherical  top  load. 


The  Q  factor  for  Spherical  Top  Load 


Figure  15.  Q-factor  for  spherical  top  load. 
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7.  CONCLUSION 


The  electromagnetic  properties  of  an  antenna  can  be  described  as  a  linear  combination  of 
eigemnodes.  The  quasi-static  antenna  design  algorithm  is  a  method  for  calculating  the  radiation 
resistance  of  the  dipole  eigenmode  (the  of  term)  and  the  DC  capacitance  of  the  antenna.  The  three- 
element  circuit  model  of  the  dipole  eigenmode  impedance  does  not  accurately  explain  the  impedance 
and  Q.  An  accurate  model  for  the  impedance  requires  a  contribution  from  the  octupole  and  higher 
eigemnodes.  The  octupole  eigenmode  can  be  approximated  with  a  parallel  capacitor.  The  interaction 
of  the  dipole  and  octupole  eigemnodes  introduces  an  of  term  in  the  radiation  resistance.  The 
octupole  eigenmode  reduces  the  antenna  Q  below  the  value  predicted  by  the  quasi-static  antenna 
design  algorithm.  The  Q  limits  based  optical  theorem  and  a  dipole  excitation  are  valid  only  in  the 
low-frequency  limit  where  the  octupole  eigenmode  does  not  contribute  to  the  radiation  resistance.. 

The  quasi-static  antenna  design  algorithm  models  a  top  load  as  a  sum  of  multipole  moments.  The 
algorithm  appears  to  converge  in  Q  and  shape.  In  the  case  of  a  spherical  enclosure,  the  Q  is  only 
slightly  larger  than  the  minimum  Q  spherical  cap  monopole.  This  algorithm  allows  antennas  to  be 
designed  with  higher  radiation  resistance  and  a  slightly  higher  Q  [7].  With  minor  modifications,  the 
same  algorithm  can  be  used  with  different  enclosing  surfaces. 
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